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ABSTRACT 

Aims. A recent research shows that particles with a spectrum of a relativistic Maxwellian plus a high-energy tail can be accelerated by relativistic 
collisionless shocks. We investigate the possibility of the high-energy particles with this new spectrum injected in pulsar wind nebulae (PWNe) 
from the terminate shock based on the study of multiwavelength emission from PWNe. 

Methods. The dynamics of a supernova remnant (SNR) and multiband nonthermal emission from the PWN inside the remnant are investigated 
using a dynamical model with electrons/positrons injected with the new spectrum. In this model, the dynamical and radiative evolution of a 
pulsar wind nebula in a non-radiative supernova remnant can be self-consistently described. 

Results. This model is applied to the three composite SNRs, G0.9+0.1, MSH 15-52, G338. 3-0.0, and the multiband observed emission from 
the three PWNe can be well reproduced. 

Conclusions. Our studies on the three remnant provide evidence for the new spectrum of the particles, which are accelerated by the terminate 
shock, injected into a PWN. 

Key words, gamma rays: theory - ISM: individual (GO.9+0.1, MSH 15-52, G338. 3-0.0) - supernova remnant 



1. Introduction 

PWNe, which are prominent sites of high-energy emission 
in the Galaxy, are powered by pulsars associated with them. 
A pulsar inside a PWN loses its rotational energy through a 
pulsar wind composed of magne t ic flux and high-energy par - 
ticles dGoldreich & Julianl 1 19691 iGelfand et all I2007L l2Q09h . 
The ultra-relativistic wind flows relativistically into a non- 
relativistic ejecta of the ambient supernova remnant (SNR), 
which results in the PWN and a termin ation shock (TS), where 
the p l asma is dec e lerate d and heated ([Reynolds & Chevalier . 
Il984 : Volpi et al. . 2008 ). High-energy particles are injected 
into the nebula from the TS, and multiband nonthermal pho- 
tons with energies ranging from radio, X-ray to 7-ray bands 
are emitted during the evolution of the PWN. 

Usually, multi-wavelength observational results of a PWN 
cannot be well reproduced by the radiation of particles in- 
jected with a single power-law spectrum, and a broken power- 
law must be employed to better explain the obse r vations 
(e.g., [A toyan & Aharonian , 1996; Aharonia n et all |2005a; 
Zha ng et all 2008 ). However, the physics behind the broken 



power-law is unclear. For a typical PWN, the wind from the 
pulsar can flow relativistically with a lorentz factor of ~ 10 6 . 
The resulting TS can accelerate particles to relativistic energy. 
On the other hand, based on the long-term two dimensional 



particle-in-cell simulations, Spitkovsky (2008) found that the 
particle spectrum downstream of a relativistic shock consists of 
two components: a relativistic Maxwellian and a high-energy 
power-law tail with an index of —2.4 ± 0.1. Based on this 
finding, with the assumption that the high-energy particles in 
a PWN are injected with a spectrum of a relativistic Maxellian 
plus a power-law high-energy tail, we investigate the multi- 
wavelength emission from PWNe to test the possibility of par- 
ticles with the new spectrum injected into the PWNe. 

The dynamical evolution of the PWN is calculat ed basi- 
cally according to the model in IGelfand et al. (2009), which 
can self-consistently describe the dynamical and radiative evo- 
lution of a pulsar wind nebula in a non-rad iative supernova 



remnant. Different from Gelfa nd et al. (2009), in which a sin- 
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gle power-law injection spectrum for the electrons/positrons is 
employed to discuss the radiative properties during different 
phase of the PWN, we argue in this paper that the high-energy 
particles are injected with the new spectrum of a relativistic 
Maxellian plus a power-law high-energy tail during the evo- 
lution in this paper, and a kinetic equation is used to obtain 
the energy distribution of the particles. We apply the model 
to the PWNe in the composite SNRs, GO.9+0.1, MSH 15-52 
and G338. 3-0.0, which have been observed in radio, X-rays and 
very high-energy (VHE) 7-rays. 

GO.9+0.1 has a 2' PWN inside a 8' shell in the radio band 



dHelfand & Becked [1987). The PWN is powered by an en- 
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ergetic pulsar PSR J 1747-2809, which was recently discov- 
ered in GO.9+0.1 with the NRAO Green Bank Telescope at 
2GHz dCamilo et all l2009l) . A jet-like feature of the nebula 
was revealed in the hi gh-resolution observation with Chandra 



(Ga ensler et al.L 120011) . The observation in the TeV band with 
HES S indicates that the P WN is a weak emitter in VHE 7-rays 
dAharonian et alll2005al) . 

MSH 15-52 is a complex SNR with a ragged shell in the 



radio observations dCaswell et all 119811) . An elongated PWN 



powered by an energetic pulsar was f ound in the remnant i n 
the X-ray o bservations with R OSAT ( Trussoni et all fl9 96). 
Bepp oSAX dMineo et all l200ll) and Chandra dGaensler et all 
20021) . The remna nt has also been well d etected in VHE 7- 
rays with HESS dAharonian et all l2005bl) and CANGAROO 



(Nak amori et all |2008), and the significant VHE 7-rays are 
identified to be produced from the PWN. 

A VHE 7-ray source HESS J1640-465 was discovered 



by HE SS in the survey of the inner Galaxy dAharonian et al 



2006), and it is spatially coincident with the composite SNR 
G338. 3-0.0. An extended X-ray source was found to be located 
at the center of th e VHE 7-ray s ource with XMM-Ne wton 
( Fun k et all 120071) . Very recently, iLemiere et al ] d20Q§ pre- 
sented the high resolution X-ray observations on the PWN, and 
a point like source as a putative pulsar appears in the X-ray 
observations. 

In this paper, we investigate the possibility of the particles 
with the new spectrum injected in PWNe based on applications 
with the spectrum to the three composite SNRs. In Section [2] 
the model is simply described, and the results from the applica- 
tions to the three SNRs GO.9+0. 1, MSH 15-52 and G338. 3-0.0, 
are presented. The main conclusions and discussion are given 
in Section[3] 

2. Model and Results 

A PWN is powered by the its pulsar which dissipates the ro- 
tational energy into the nebula. The spin-down luminosity of a 
neutr on star with a rotation period of P ev olves with time as 
(e.g. JGaensler & SlaneL l2006: Slanel l2008l) 



E(t) =E Q (l + t/ TQ ) 



-(n+l)/(n-l) 



(1) 



where ro is the spin-down time scale of the star, Eq is the initial 
spin-down power, n is the braking index of the pulsar, which is 
equal to 3 for magnetic dipole spin-down. 

Magnetic field and high-energy particles are injected at 
the TS located where the ram pressure of the unshocked 
wind is equal to that of the PWN. In this paper, we as- 
sume the spin-down power is distributed between electrons 
and positions (E e (t) = r] e E(t)), a nd magnetic fields (Eb = 
r) B E(t)) (e.g.. IGelfand et all l2009h . In IGelfand et"aD (I2009I) . 
they used a single power-law injection spectrum for the elec- 
trons/positrons to discuss the radiative properties during dif- 
ferent phases of the PWN evolution. However, a broken 
power-law spectrum is usually needed to reproduce the non- 
therm al emission of a PWN wit h multi-band observations 
(e.g., Atovan & Aharonian, 1996; Venter & de Jager, 2006; 
Slane et all 120081; IZhang et all 120081) . Recently, based on 



Smtkovsky (2008) found that the particle spectrum down- 
stream of a relativistic shock can be fitted as a Maxwellian plus 
a power-law tail. The form is 7(7) = C17 exp(— 7/A71) + 
C27 _a min[l, exp(— 7 — 7 C ut)/A7 C ut], where 7 is the Lorentz 
factor of the particles, A71 is close to (70 — l)/2, 70 is the 
Lorentz factor of the upstream flo w, Ci = for 7 < 7 m i n , 
7 m ; n is about seven times of A71 (Snitkovskvl l2008h . In this 
paper, we investigate the possibility of particles with this spec- 
trum injected in PWNe by multiwavelength studies. 

We argued that high-energy particles injected in a PWN are 
accelerated by the TS, for which the upstream pulsar wind typ- 
ically can have a Lorentz factor ~ 10 6 . Therefore, we assume 
the spectrum of the high-energy particles injected in the PWN 
has the form, 



Q(E,t) 



C(*)#exp 



C(t) 



(-*) 



E < E„ 



if ex P 



-&)+/(& 



<E < E n 



where, a = 2.4 ± 0.1, E h ~ 2.6 x 10 5 7 ts , 6 MeV, 7 ts , 6 is 
the Lorentz factor of the upstream pulsar wind of the TS in 
units of 10 6 , E min = f Blin E h with / min ~ 7, / is normal- 
ized by E m 3 a /E h exp(-E min /E b ) j = f (E m -, T1 /E h )~ a . This 



spectrum is similar to that given in Snit kovskvl 0008) except 
for a discontinuity at E m - ln , which has no influence on the final 
multiband photon spectrum, and the cut-off form at the high- 
est energy. We leave E max as a parameter constrained by the 
observational results for a PWN. C(t) can be obtained with 



C{t) 



E e (t) 



















(V) 2 1 



(3) 



Assuming the particles are homogeneous distributed in space in 
the PWN, and the distribution in the emission region is given 
by 



dN(E,t) 
di 



d_ 

dE 



EN(E, t) +Q(E,t), 



(4) 



where E is the energy-loss rate of the particles with energy E. 
Energy-loss mechanisms include synchrotron radiation, inverse 
Compton scattering and adiabatic losses. 

The dynamics of the PWN in the supernova shell is calcu 



the long-term two-dimensional particle-in-cell simulations, 



lated b asically following the model presented in Gelfa nd et al 
d2009h . The model assumes that the progenitor supernova 
ejects material with mass M C j and energy E sn into the ambi- 
ent matter with a constant density pism- Assuming the PWN 
has no influence on the dynamics of the forward shock and the 
reverse shock, the velocity and the radius of the forward shock 
and the reverse shock of the surrounding SNR a re calculated 
with the equations in iTruelove & M cKee dl999l) . The pulsar 
dissipates energy into the PWN, which sweeps up the super- 
nova ejecta into a thin shell surrounding the nebula, and new 
particles are injected into the PWN at each time step. 

We now apply the model the investigate the three compos- 
ite SNRs, GO.9+0.1, MSH 15-52 and G338. 3-0.0. Firstly, if the 
distance, the age and the radii of the SNR shell and the PWN 
are known, the parameters such as the supernova energy (E an ), 
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Fig. 1. Upper panel: radius of PWN (i? pwn ) for M e j = 2 
(dashed line), M e j = 6 (dotted line), M e j = 8 (solid line) and 
M e j = 10 (dash-dotted line) Mq, respectively, with the other 
parameters listed in Table Q] for GO.9+0.1. Lower panel: mag- 
netic field strength in the PWN (B pwn ) with different M e j, and 
the others are same as the upper panel. 



the ejecta mass (M e j), the ambient density (pismX an d the ini- 
tial spin-down energy (Eq) can be constrained by making the 
results consistent with the known values. In our calculation, 
a typical spin-down time scale (r S d) is set to 500 yr for the 
three remnants, and Eq is also constrained by the current spin- 
down power of the pulsar if it is available. Finally, in order to 
obtain a consistent nonthermal emmission with the multiband 
observational fluxes for the SNR, the other parameters can be 
constrained. Note that uncertainties for some parameters still 
exist due to both their insensitivities on the final results and the 
uncertainties of the key properties of the system, such as the 
distance, the age and the properties of the pulsar. 

2.1. GO.9+0.1 

The composite SNR GO.9+0.1 con sists of a bright 2' PWN 
hosted by a 8' shell in the radio band dHelfand & BeckerMl987 ; 



La Rosa et all 2000). Recent high resolution radio study of 
the PWN indicates that the radio emission of the nebula has 
a spectral index of —0.18 ± 0.04, and the fluxes are 1.35, 1.45 
and 1 .75 Jy for the w avelength 3.6, 6 and 20 cm, respectively 



( Dub ner et al.L 120081) . Furthermore, an energetic pulsar PSR 



in GO.9+0.1 w ith the NRAO Green Bank Telescope at 2GHz 
dCamiloetalll2009t) . 

The PWN in GO.9+0.1 was detected in X-rays with 
BeppoSAX, and the result shows that the X-ray emission from 
the source has a power-law spectrum with a photon i ndex of 
2.0 ± 0.3 dMereghetti et all 11998k ISidoli et all l2000h . High- 
resolution observation with Chandra indicated that the neb- 
ula is axisymmetric with a jet-like feature, which is evidence 
of a torus of emission in the pulsar's equatorial plane and a 



jet directed a l ong th e pulsar spin axis ( Ga ensler et al. , l200ll) . 
IProquet et all (2003) studied GO.9+0.1 using observations by 
XMM-Newton. The X-ray spectrum softens with distance from 
the core, and the spectrum in the energy range 2-10 keV has 
a power law form with a photon index of ~ 1.9. 

Very High-energ y (VFJE) emission from G O.9+0. 1 has been 
detected with HESS dAharonian et alll2005al) . The photon flux 



above 0.2 TeV is 5.7 x 10~ 12 cm~ 2 s _1 , and the spectrum can 
be fitted with a power-law with a photon index 2.4 ± 0.31. The 
source is a weak TeV emitter, and the VHE 7-rays appear to 



origina te from the core rather than the shell dAharonian et al 
l2005ah . 

The dynamical and radiative properties of the composite 
SNR GO.9+0.1 are investigated with the parameters in Table 
□ for this source. Although t he pulsar PSR J 1747 -2809 has a 
characteristic age of 5.3 kvr. ICamilo et al. I d2009h argued that 
GO.9+0.1 has a small age of no more than 2 — 3 kyr either 



from PWN evolution models of iBlondin. Chevalier & Frierson 
(2001) for the o bserved ratio R uvin / R sn i- — 0.25 or from the 
PWN energetics dDubner et alll200 8). The distance of the pul- 
sar is likely in the range of 8 kpc to 16 kpc due to the uncer- 
tainty of the ele ctron density model toward the distant inner 



Galactic regions dDubner et a l.. 2008). We assume the distance 



is 8.5 kpc in the calculation, then the radii of the PWN and 
the shell are 2.55 pc and 10.2 pc, respectively. Moreover, with 
E sn = 10 51 erg and M ej = 8M , an age of 1900 yr and a 
relatively low density n- ism = 0.01 cm~ 3 are needed to well 
reproduce structure of the system, i.e., the radius of the shell 
and the ratio i? pwn / i? S nr- A pulsar's velo city of 120 km/s sim 



J1747-2809 with a period of 52 ms has recently discovered 



ilar as the Crab pulsar (Kapl an et al. . 2008 ) is used to illustrate 
the evolution of the PWN, which has no influence on the final 
results for GO.9+0.1 since it is a young remnant, and now the 
pulsar is safely in the nebula. With these parameters, the re- 
sulting radii of the PWN and the SNR shell are 2.6 and 10.2 
pc, respectively, consistent with the observational results (up- 
per panel FigJTJ. The influence of the ejecta mass M G j on the 
radius of the nebula is shown in the upper panel of FigQ] With 
a smaller M e j, both the nebula and the SNR shell expand more 
quickly, and the nebula collides with the reverse shock at later 
time. As a result, the magnetic field in the PWN is weaker for a 
smaller M e j due to the bigger volume of the nebula (the lower 
panel in FigJTJ. 

The dynamical properties of the SNR during the evolution 
is similar as those for the Crab remnant in Gelfand et al. I d2009h 
although the new spectrum of the particles is used in this paper. 
Initially, the pressure of the PWN is much bigger that of the 
surrounding supernova ejecta, so the PWN expands adiabat- 
ically into the cold supernova ejecta. The ejecta surrounding 
the PWN is swept up to a thin shell, which is decelerated by 
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Table 1. Input Parameters for the three PWNe 



Parameter 


GO. 9+0.1 


MSH 15-52 


G338. 3-0.0 (CASE 1) 


G338. 3-0.0 (CASE 2) 


Distance (kpc) 


8.5 


5.2 


10.0 


10.0 


T ago (yr) 


1900 


1900 


4500 


8200 


Tsd (yr) 


500 


500 


500 


500 


M ej (M ) 


8.0 


8.0 


8.0 


8.0 


E sn (10 51 erg) 


1.0 


3.0 


0.5 


0.8 


n ism (cm -3 ) 


0.01 


0.01 


0.1 


1.0 


E (10 39 erg s" 1 ) 


2.0 


1.5 


10 


0.5 


n 


3.0 


2.84 


3.0 


3.0 


m (io -3 ) 


1.0 


20.0 


1.0 


0.3 


Q 


2.5 


2.4 


2.5 


2.5 


S max (TeV) 


900 


1000 


500 


500 


E h (10 5 MeV) 


0.2 


0.9 


0.8 


0.8 




Electron Energy (MeV) Photon Energy (MeV) 



Fig. 2. Left panels: particle spectra of electrons and positrons at different times during the evolution of the PWN. Right panels: 
the resulting spectral energy distribution at different ages with the other parameters listed in TableQ] 



ram pr essure since its velo city is higher than the local sound 
speed dGelfand et all 120071) . The mass of the PWN M pwn in- 
creases continuously since the shell expands faster than the am- 
bient ejecta. This expansion phase ends when the PWN col- 
lides with the reverse shock of the SNR. After the collision, 
the pressure of the nebula P pwn is much less than that of the 
supernova ejecta around the nebula P S nr(^?pwn)- The velocity 
of the PWN decreases greatly, and finally the PWN is com- 
pressed when the shell moves inwardly. During this process 
of compression, the magnetic field strength in the nebula in- 



creases significantly (lower panel in FigJTJ, and, as a result, the 
synchrotron luminosity has a rapid rise (lower right panel in 
Fig |2]i . Furthermore, the radius of the PWN decreases signifi- 
cantly, and the PWN will expand again when the pressure in it 
eventually becomes bigger than that of the surrounding ejecta. 
The nebula enjoys a series of contractions and re-expansions 
until the SNR enters the radiative phase of its evolution. The 
pulsar with a velocity moving in the space will leave the PWN 
if the velocity is high enough in the compression process, and 
it can re-enter the nebula when the nebula expands again. 
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Particles with a spectrum of a Maxwellian plus a power-law 
tail (a = 2.5, / = 1.0) are injected in the PWN when the pul- 
sar is inside it. The particles lose energy through synchrotron 
radiation, inverse Compton scattering and adiabatic loss. Non- 
thermal emission from the PWN from the radio to X-ray band 
is produced via synchrotron radiation, whereas 7-rays are pro- 
duced through inverse Compton scattering on the seed photons, 
i.e., cosmic microwave background (CMB), infrared (IR) and 
optical (opt) photons. For the soft seed photons, the energy den- 
sity of the infrared radiation is 0.23 eV cm" 3 as used fo r the 
GALPROP code dStrong et all bOOOt IPorter et all 120061) and 
a density of 5.7 eV cm~ 3 for the optical component, which 
is 50% smaller than the value used in GALPROP, are used in 



Aharonian et al.l d2005al) to discuss the origin of the VHE 7- 
rays. We find densities of 0.5 eV cm -3 for the infrared compo- 
nent (fcTiR = 3 x 10~ 3 eV) and 20 eV cm~ 3 for the optical 
soft photons (fcT op t = 0.25 eV) can reproduce the observed 
spectrum in the VHE 7-ray band well, so these densities are 
used in this paper. The resulting energy distribution of the par- 
ticles in the PWN and the multi-wavelength non-thermal emis- 
sion during the evolution are presented in Figj2] The particles 
with energies smaller than Et, are mainly from the high-energy 
particles with energies > Et, experienced energy losses, and 
the energy distributions in the range < Et, at different age in 
the evolution all show a power-law form with an index of ~ 1 
(dN/dE oc E~ 1 ). On the other hand, the energy distribution in 
the higher-energy band (> Et,) becomes steady between 1000 
yr and 20000 yr, and the distribution can also be represented as 
a power-law but with an index ~ 2.5. 

The expansion and compression of the PWN cause the 
magnetic field strength in the nebula B pwn to decrease and 
increase, respectively. The synchrotron emissivity diminishes 
before 10000 yr as the PWN expands into the cold supernova 
ejecta, when the magnetic field strength in the nebula decreases 
gradually. However, synchrotron radiation can be important 
again when the PWN is compressed, because the magnetic field 
strength rises greatly during the compression. In this case, the 
energy of the particles in the nebula can have a rise due to the 
adiabatic compression, then the emissivity of inverse Compton 
scattering can also have a rise (see the emission in Fig|2]for 
30000 yr). During the whole evolution process, the PWN is 
an important emitter of the GeV 7-rays, although the emission 
from the radio to the X-ray band is insignificant sometimes. 

Fig[3] shows the evolution of the energy contained in the 
PWN (E pwn ) during the evolution of the system, and the spin- 
down power of the pulsar and the adiabatic loss rate of the 



PWN, which can be obtained with L 



ad 



EL 



are indicated in Fig|4] Here v pwn and R pwn are the velocity 
and radius of the PWN, respectively. At first, the energy of 
PWN increases continuously since the pulsar injects the spin- 
down power into the nebula, and this stage ends after ~ 1000 
yr (Figj3j) when the adiabatic loss rate is considerable com- 
pared with the spin-down power of the pulsar (Fig0]i. After 
this phase, the energy of the nebula decreases or increases due 
to the competition between the spin-down power of the pulsar 
adiabatic loss. 

With the parameters in Table[T] the dynamical structure and 
the multi-band radiative properties of GO. 9+0. 1 are well repro- 




10 1 



10 10 10 

Time (Years) 



10 s 



Fig. 3. The energy contained in the PWN at different ages of 
the system with the parameters same as Fig|2] 




10 10 10 

Time (Years) 



10 



Fig. 4. The spin-down power of the pulsar (solid line) and the 
adiabatic loss rate of the PWN (dotted line) during the evolu- 
tion of the system. The others are the same as Figj3] 



duced at 1900 yr. At this age, the magnetic field strength in the 
PWN is 8. 1 fiG, and the current spin-down power of the pulsar 
is 8.67 x 10 37 erg s _1 , whi ch is still consistent with the mea- 
sured one 4.3 x 10 37 erg s _1 ( Camilo et al. , 20091) due to the un- 
certainty of the moment of inertia; moreover, the observational 
radii of the nebula and the SNR shell are reproduced in this sce- 
nario. The resulting multi-band emission of the PWN from the 
model is shown in the lower panel in Figj5] The radio and the 
X-ray observations can be well explained as synchrotron radia- 
tion of the particles injected in the PWN; VHE 7-rays from the 
nebula are mainly produced by inverse Compton scattering on 
the optical light and the cosmic microwave background. 

With a bigger ejecta mass (M j), the radius of the PWN is 
smaller due to the deceleration of the ejecta matter before the 
collision between the reverse shock and the nebula. As a result, 
the resulting synchrotron radiation is stronger for a bigger M e j 
(the upper panel in Fig|5]). Moreover, the influence of the break 
energy Et, on the resulting multiband nonthermal emission is 
indicated in Fig|6]for E h = 0.1 x 10 5 (dashed line), 0.2 x 10 5 
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Fig. 5. Upper panel: the multiband spectral distribution of the 
nonthermal emission for M e j = 2 (dashed line), M e j = 6 (dot- 
ted line), M C j = 8 (solid line) and M e j = 10 (dash-dotted line) 
Mq with the other parameters listed in Table Q] for GO.9+0.1. 
Lower panel: the resulting photon emission of synchrotron 
(dashed line), inverse Compton scattering on the CMB (dot- 
ted line), IR(dash-dotted line), starlight (short-dashed line), and 
synchrotron photons (close-dotted line) for M e j = 8M and 
the others same as the upper panel. The solid line represents 
the whole inverse Compton scatteri ng of the electrons/ positions 
on all the soft photon s. The radio jDubner et all I2008J). X-ray 
dProquet et al.ll2003l) and VHE ^-ravflAharonia n et all [2p05a) 
observations for the PWN are also shown in the figure for com- 
parison. 




10 10 10 10 10 10 10" 10 10 10 10 10 

Photon Energy (MeV) 

Fig. 6. The resulting spectral energy distributions for E\, = 
0.1 x 10 5 (dashed line), 0.2 x 10 5 (solid line), 0.5 x 10 5 (dotted 
line) and 1.0 x 10 5 (dash-dotted line) MeV, respectively. The 
others are the same as the lower panel in Fig[5] 




10" 13 10" 11 10" 9 10" 7 10" 5 10" 3 10"' 10' 10 3 10 5 10 7 10 9 

Photon Energy (MeV) 

Fig. 7. Comparison of the resulting photon emission of syn- 
chrotron (dashed line), inverse Compton scattering on the CMB 
(dotted line), IR(dash-dotted line), starlight (short-dashed line), 
and synchrotron phot ons (clo se-dotted line) with the radio 
(IGaensler etall l2002h X-rav jMineo et all l200ll) and VHE 
T-rav dAharonian et a l., 2005b) observations for the PWN in 
MSH 15-52. The parameters are listed in TabkQ]for this rem- 
nant. 



(solid line), 0.5 x 10 5 (dotted line) and 1.0 x 10 5 (dash-dotted 
line) MeV, respectively. With a bigger E\>, the resulting fluxes 
in X-rays and 7-rays are higher. The multiband observational 
results in the radio, X-rays and VHE 7-rays can be well re- 
produced with _Eb = 0.2 x 10 5 . In such a case, the Lorentz 
factor of the pulsar wind upstream of the TS can be estimated 
as - 0.8 x 10 5 . 

2.2. MSH 15-52 

MSH 15-52 (G320.4-1.2) is a complex object located at a dis- 
tance of 5.2 ± 1.4 kpc based on the HI absorption measurement, 
consistent with th e value of 5.9 ± 0.6 kpc from the pulsar dis- 
persion measure (ITavlor & Cordesl 119931) . It is a rough circu- 



lar SNR with a diameter of ~ 30' from the radio observations 



(ICaswell et all ll981). and a PWN powered by an energetic pul 



sar PSR B1509 -58 (e.g.. IGaensler et alll2002b was discovered 
in the remnant. iLivingstone et al. I d2005l) presented an updated 
timing solution for the young energetic pulsar PSR B 1509-58 
based on the 21.3 yr of radio data and the 7.6 yr of X-ray tim- 
ing data. The results show that the pulsar has a period of about 
150 ms with a braking index of n — 2.839 ± 0.003, and a 
characteristic age of ~ 1700 yr. In the radio band, the source 
appears as a shell-lik e SNR with bright spots in the NW and SE 
region (IMineo et all 1200 11) . and the N W region co i ncide s with 



the IR and optical nebula RSW 89 (ISeward et all [1982). The 



PW N around the pulsar i s comparatively faint in radio bands, 
and IGaensler et all d2002l) estimated a flux density of 2 ± 1 Jy 
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at both 0.8 and 1.4 GHz u s ing th e data of IWhiteoak & Greenl 
19961) and iGaensler et aD {[999) for the PWN. In the X-ray 



band, diffuse emission is detected from the shell-like remnant 
with a bright spot coincident with the NW zone. The PWN has 
an elongated structure roughly centered on the pulsar with two 
arms extending several arcminutes along the N W and SE direc 



tions in the X-ray ob servations with ROSA T (iTrussoni et al 



1996b and Chandra dGaensler et all l2002h . The hard X-ray 
spectrum for the PWN was measured by BeppoSAX and the 
power-law fit of the unpulsed high energy flux gave a photon in- 
dex 2.1 ±0.2 and a flux ofJJ2>£l£-i°erg cm " 2 s" 1 in the 20- 



200 keV energy range (Mineo et al., 2001). At VHE 7-ray en- 



ergies, the SNR MSH 15 -52 has been observe d with the H ESS 



( Ahar onian et all l2005bl) and CANGAROO dNakamori et al 



2008|). The TeV 7-ray s are from an elliptically shaped region, 
and the jet extends more prominently to the south/southeast. 
This morphology co incides with the diffuse PWN as observed 
at the X-ray bands ( Aharonian et al., 2005b]). The overall en- 
ergy spectrum can be fitted by a power-law with photon index 
T = 2.27 ± 0.03 sta t ± 0.02 syst in the energy range 280 GeV 
to 40 TeV in the HESS result, and a compatible result has been 
obtained with the CANGAROO telescope. 




Photon Energy (MeV) 

Fig. 8. Comparison of the resulting photon emission of syn- 
chrotron (dashed line), inverse Compton scattering on the CMB 
(dotted line), IR(dash-dotted line), starlight (short-dashed line), 
and synchrotron pho tons (close-dotte d line) with the X-ra 
dLemiereet aUl2009t ) and VHE 7-ray dAharonian et alll200 



<6) 



observations for the PWN in G338. 3-0.0. The parameters are 
listed in TabliQfor G338.3-0.0 (CASE 1). 



The age of the SNR MSH 15-52 has been estimated as 6- 
20 kyr with the standard paramete rs for the inter s tellar medium 
and for the supernova explosion dSeward et all 1 1 9 8 3D . which 



is significantly larger than the characteristic age (1700 yr) of 
pulsar. However, it is very likely that the system is young with 
an age of ~ 1700 yr since the SNR has e xpanded rapidly into 



a low-density cavity dSeward et al.L 119831) . whic h is supported 



by the observation of HI emission in this region dDubner et al 
20021) . 



The resulting multiband nonthermal emission from MSH 
15-52 are shown in Fig|7l and the parameters are listed in 
Tabl^U With an initial spin-down power of 1.5 x 10 39 erg 
s _1 and To = 500 yr, the spin-down power of the pulsar is 
5.7 x 10 37 erg s _1 at an age of 1900 yr, consistent with the 
measured one, 5.4 x 10 37 erg s -1 dLivingstone et alll2005i) . An 
energy density of 1.5 eV cm~ 3 for the optical soft photons is 
used in the calculation , which is consisten t with the value from 
the GALPROP code dStrong et all l2000t IPorter et all 120061) . 
For the infrared component, we find a density of 1.5 eV cm~ 3 
is needed to well reproduce the observation in VHE 7-rays, so 
this value is employed in the calculation. With the parameters 
in TabldTJfor MSH 15-52, the radii of the PWN and the SNR 
shell are 3.6 pc and 14.8 pc, respectively, and the resulting mag- 
netic field strength inside the PWN is 19.3 uG now . The ob- 
served emission from the radio (Gaensler et al., 2002) to the X- 
ray bands detected with BeppoSAX dMineo et alluOOll) can be 
well explained as the synchrotron radiation of the high-energy 
electrons/positrons injected in the nebula; moreover, the result- 
ing emission from 100 MeV to 10 TeV is mainly produced via 
inverse Compton scattering off the soft infrared photons, and 
the HESS flux points can be well reproduced. 



2.3. G338. 3-0.0 

HESS J1640-465 was discovered in the survey of the inner 
Galaxy with HESS as a center filled VHE 7-ray source with 
a differential energy spectrum fitted as a power-law with a n 
index of 2.42 ± 0.14 above 0.2 TeV dAharonian et all 120061) . 
This source is spatially coincident with the SNR G338.3- 
0.0, which has a broken shell with a diameter of 8' indi- 
cated in the the 843 MHz radio survey using the Molongl o 
Observatory Synthesis Telescope ( Whiteoak & Greenl 19961) . 
A highly absorbed source inside t he shell was indicate d in 



the X- ray observat i ons w ith ASCA (ISugizaki et all 12001) and 



Swift ( Landi et al. , 2006). Extended X-ray emission centered 
about the VHE 7-rays was dete cted with a dedicated XMM- 
Newton observation in 2006 ( Funk et al. . 20071) . Recently, 



Lemie re et al. presented the high resolution X-ray ob- 

servations with Chandra on the extended source. The observed 
morphology shows a PWN and a point-source presented as a 
potential pulsar. The spectrum of the putative pulsar and the 
PWN can be fitted with a power - law w ith an index of 1 . 1 and 



2.5, respectively (ILemiere et al.L 120091) . They argued that the 
pulsar's spin power is ~ 4 x 10 36 erg s _1 based on the X- 
ray luminosities of the putative pulsar and nebula between 2 to 
10 keV. The distance of the source is less certain, and it can 
be from 8 kpc to 13 kpc based o n the H I absorption features 



observed along the line of sight (ILemiere et al., 2009) 



Assuming the SNR has a dis tance of 10 kpc, the radi us of 
the radio shell is - 12 pc (4') dWhiteoak & Greenl fl996l) . and 
that of the PWN is ~ 3.5 pc (1.2') in the X-ra ys according to 



the hi gh-resolution observations with Chandra (ILemiere et al 
2009). Firstly, G338.3-0.0 is assumed to be a young SNR ex- 



panding into a tenuous medium with a density of 0.1 cm~ 3 . 
The age of the remnant should not be too young since the PWN 
must contain enough energetic particles to produce the obser- 
vational VHE 7-rays through inverse Compton scattering. As 
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Photon Energy (MeV) 



Fig. 9. The resulting spectral energy distribution with the pa- 
rameters listed in TabldT]for the CASE 2 of G338. 3-0.0 (CASE 
2). Others are the same as FigJS] 



a result, the supernova explosion energy E sn and the initial 
spin-down power of the pulsar are chosen to be 10 51 erg and 
10 40 erg s _1 , respectively; moreover, E max is constrained to 
500 TeV to make the resulting emission consist with the ob- 
servations in X-rays, and relatively high densities are needed 
to reproduce the VHE 7-ray fluxes detected with H.E.S.S., i.e., 
5.0 eV cm -3 and 25.0 eV cm~ 3 for the infrared and the optical 
soft photons, respectively. At an age of 4500 yr, with the other 
parameters listed in TabldT]for the SNR G338.3-0.0 (CASE 1), 
the radii of the shell, the reverse shock and the PWN are 6.9, 8.8 
and 1 1 . 1 pc, respectively, and the magnetic field strength in the 
nebula is 3.7 /iG. In such a case, the pulsar is energetic with a 
spin-down power of 1.0 x 10 38 erg s" 1 now. F urthermore, the 



obser ved spectra in the X-rays with Chandr a dLemiere et al 



2009) and in the VHE 7-rays with H.E.S.S. (Ah aronian et al 



2006) can be reproduced (Fig|HJ. Inverse Compton scattering 



on the optical light dominates the resulting emission at 0. 1 TeV, 
whereas the emission above 1 TeV is mainly produced via in- 
verse Compton scattering off the IR photons. 

On the other hand, we investigate the dynamical and radia- 
tive properties of G338. 3-0.0 expanding into a medium with a 
density of 1 cm~ 3 . With E c - } = 0.8 x 10 51 erg, the radius of the 
shell of the SNR in radio means that the age of the remnant is 
~ 8000 yr. With an initial spin-down power of 0.3 x 10 39 erg 
s _1 and other appropriate parameters (see Tabla[T]for G338.3- 
0.0 CASE 2), the current spin-down power of the pulsar is 
1.65 x 10 37 erg s _1 , and the resulting radius of the PWN is 
just 2.45 pc, which is even smaller than the exte nsion of the X 
rays ( ~ 3.5 pc) assuming a distance of 10 kpc (ILemiere et al 
120091) . A smaller rye = 0.3 x 10~ 3 , and the maximum energy 
of the particles is set to 500 TeV to reproduce the observational 
fluxes in the X-rays and 7-rays. The resulting multiband non- 
thermal emission is indicated in Fig(9]with soft densities of 1.0 
eV cm~ 3 and 6.0 eV cm~ 3 for the infrared and the optical soft 
photons, respectively. In this scenario, the PWN has been com- 
pressed by the reverse shock, and the resulting flux with ener- 
gies below 1 eV is about two orders of magnitude higher than 
that in the CASE 1. 



Motivated by the finding that the spectrum of the particles 
downstream of a relativistic shock consists of two components: 
a relativistic Maxwellian a nd a power-law hig h-energy tail with 
an index of —2.4 ± 0.1 ( Spitk ovskvt 120081) . we investigate 
the possibility of particles with this new spectrum injected in 
PWNe from the TS based on the studies of multiband emis- 
sion from PWNe. F ollowing the dynamical method proposed in 
Gelfan d et al. (2009), we study the dynamical and multi-band 



radiative properties of the three composite SNRs GO.9+0.1, 
MSH 15-52 and G338. 3-0.0. With appropriate parameters, we 
find that a typical PWN is an important 7-ray emitter during 
its evolution although the non-thermal radiation from the ra- 
dio to the X-ray band is insignificant sometimes. The multi- 
band observations of the three PWNe in the remnants can be 
well reproduced with the new spectrum of the injected parti- 
cles. Therefore, our studies on the dynamical and multiwave- 
length radiative properties of PWNe provide evidence of high- 
energy electrons/positrons can be injected into a PWN with a 
Maxwellian plus a power-law high-energy tail from the TS of 
the PWN. 

In modeling the multiband nonthermal emission from a 
PWN detected in the radio, X-ray and 7-ray bands, particles 
injected with a spectrum of a broken power-law are widely 
used to reproduce the observed multiwavelength emission ( e.g., 
Venter & de Jageil l2006t Island, l2008t IZhang et all 120081) , Of 



course, for the three PWNe discussed in this paper, the multi- 
band observed spectra of them can also be explained if the par- 
ticles are injected with a broken power-law. However, it is un- 
clear why the broken power-law spectrum is valid when using it 
to reproduce the multiwavlength emission from a PWN. From 
our calculations, we have found out that the energy distribu- 
tion of the electrons/positrons in the nebula can be approxi- 
mated as a broken power-law with an index ~ 1 in the lower- 
energy band and an index of ~ 2.5 in the higher-energy part be- 
fore the PWN undergos significant compression, which is most 
likely the physical explanation of the broad usage of a broken 
power law in modeling the multi-band non-thermal emission 
from PWNe. 

In this paper, high-energy electrons/positrons are injected 
into the PWN from the TS, and the main energy of the nebula 
is contained in these particles. These particles undergo radia- 
tive and adiabatic losses when the nebula evolves in the host 
SNR. Our study indicates that, for a typical PWN with the pa- 
rameters similar as GO.9+0.1, the adiabatic loss of the parti- 
cles in the nebula is significant after an age of ~ 1000 yr (see 
Fig|3]and Fig|4]i. Multiwaveband nonthermal emission from a 
PWN has been investigated using a simplified time-dependent 
injection model, in which hi gh-energy electrons/positrons are 
injected into the PWN ( e.g., IVenter & de J ager. 2006; Slane, 



2008; Zh ang etal .1120081) . The pulsar inside the PWN transfers 



a part of its spin-down power to the particles with a spectrum of 
a broken power-law. In the si mplified time-dependent injection 
model in lZhang et al. (2008), synchrotron loss of the particles 
is taken into account, whereas the adiabatic one is ignored. As 
a result, either a relatively smaller initial spin-down power of 
the pulsar or a smaller efficiency of the power to the kinetic 
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energy of the accelerated electr ons/positrons is em ployed in 
the model. Moreover, note that in Zhan g et al . (2008), an initial 



spin-down power of 1 x 10 erg s for MSH 15-52 was used 
to investigate the multiband emission from the PWN, which is 
a factor of 15 smaller than that used in this paper. Besides the 
above reasons, the another main one is a relatively big spin- 
down time scale of ~ 5000 yr, which is oc Eq in the paper, 



used by Zhan g et al.1 (120081) . whereas in this paper it is adopted 
to be 500 yr. The energy released by the pulsar is mainly deter- 
mined by Eq min{r age , t q |, and the value in this paper is not 



much bigger than that in IZhang et alj (120081) . Therefore, the 



multiband observed spectra for MSH 15-52 can be reproduced 
within the two scenarios even the initial spin-down power of 
the pulsar is significantly different. 
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